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Abstract 

In this paper we present the analytic form of the heavy-quark coeffi- 
cient functions for deep-inelastic lepton-hadron scattering in the kinematical 
regime ^ . Here and stand for the masses squared of the 
virtual photon and heavy quark respectively . The calculations have been 
performed up to next-to-leading order in the strong coupling constant a<j 
using operator product expansion techniques. Apart from a check on earlier 
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calculations, which however are only accessible via large computer programs, 
the asymptotic forms of the coefficient functions are useful for charm pro- 
duction at HERA when the condition » ml is satisfied. Furthermore the 
analytical expressions can also be used when one applies the variable heavy 
fiavour scheme up to next-to-leading order in ag. 
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1 Introduction 



The study of deep inelastic electroproduction has led to important informa- 
tion on the structure of the proton. This information is extracted from the 
structure functions F2{x,Q'^) and Fl^x^Q"^) which appear in the cross sec- 
tion for deep inelastic lepton-hadron scattering. Here x denotes the Bjorken 
scaling variable and is the mass squared of the virtual photon exchanged 
between the lepton and the hadron. In the framework of perturbative Quan- 
tum Chromo dynamics (QCD) these structure functions can be described by a 
convolution of parton densities with coefficient functions. The latter are cal- 
culable order by order in perturbation theory. However the parton densities 
cannot be computed yet since they are of a nonperturbative origin and have 
to be extracted from the data. Starting from the late sixties these densities 
have been obtained from many experiments. Until recently the analysis was 
carried out in the kinematical range 0.01 < x < 0.95 and < 300 (GeV/c)^. 
However since the advent of the HERA accelerator the kinematical region 
has been extended to much smaller values of x {x > 10^^) and much larger 
values of {Q^ < 2 x 10^ (GeV/c)^). The most recent results come from 
the HI and ZEUS experiments at HERA, see and respectively. 

The low-x region is of great experimental as well as theoretical interest. 
The structure function F2{x, Q"^) rises very steeply when a; ^ which can be 
mainly attributed to a corresponding increase in the gluon density. There- 
fore this density is a very important issue in the investigation of the small-x 
structure of the proton. Since the gluon density appears together with the 
other parton densities in most cross sections, one has to look for those specific 
reactions in which it plays a dominant role so that it can be isolated from 
its partners. One of these processes is extrinsic charm production in deep in- 
elastic electron-proton scattering. Here the dominant production mechanism 
is represented by the photon-gluon fusion process which is indeed the only 
one in the Born approximation. Next-to-leading order (NLO) calculations 
[0], to which also other processes contribute, reveal that this picture remains 
unaltered. 

Apart from the interest in the gluon density, charm production also re- 
vived the important issue of how to treat the charm quark in deep inelastic 
scattering processes. Here one can distinguish between intrinsic IQ and ex- 
trinsic charm production. In the former case the charm is considered to be 
a part of the hadronic wave function and it is described by a parton density 
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in the hadron like the other hght flavours {u,d,s) and the gluon (g). How- 
ever this prescription is only correct if the charm quark can be treated as 
a massless particle which is certainly not the case at small-Q^ values where 
threshold effects become important. In this region the charm quark has to be 
treated as a massive particle. On the other hand when ^ ,where is 
the charm mass, extrinsic charm production via the photon-gluon fusion pro- 
cess and its higher order QCD corrections reveal large logarithms of the type 
ttg ln'^(Q^/m^). These large corrections bedevil the perturbation series and 
have therefore to be resummed via the renormalization group equations. This 
resummation entails the definition of a charm parton density in the hadron. 
Although a recent investigation ^ shows that the above logarithms lead to 
a rather stable cross section for charm production with respect to variations 
in the factorization and renormalization scale, the size of these large correc- 
tions warrants a special treatment. This is provided by the so-called variable 
flavour scheme (VFS) 0, [0- In this scheme the treatment of the charm 
depends on the values chosen for Q^. At low Q^- values the deep inelastic 
structure functions are described by the light parton densities {u,d,s,g). The 
charm contribution is given by the photon-gluon fusion process and its higher 
order QCD corrections. At large the charm is treated in the same way 
as the other light quarks and it is represented by a charm parton density in 
the hadron, which evolves in Q^. In the intermediate region one has to 
make a smooth connection between the two different prescriptions. In [||, 
and |jlO| this was done by adding and subtracting certain mass factorization 
terms. Notice that the above considerations also apply to bottom production 
when gets extremely large. 

Up to now the VFS-scheme has only been applied to heavy flavour electro- 
production using the Born approximation to the photon-gluon fusion process 
[|] , and [|I^ . It is our aim to extend this scheme to next-to-leading order 
in as- For that purpose we will calculate the full two- loop operator matrix 
elements containing one heavy quark loop. This calculation provides us with 
the terms containing the logarithms In*"' (/i^/m^) which have to be subtracted 
from the charm cross section so that the flnal result becomes independent of 
the charm mass rric. Here /x denotes the operator renormalization scale which 
can be identified with the factorization scale. Using the two-loop operator 
matrix elements up to non-logarithmic terms and the NLO light parton coef- 



ficient functions in ||TT[, we can construct an analytic form of the NLO heavy 
quark coefficient functions in the limit ^ m^. We will refer to these re- 
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suits as the asymptotic heavy quark coefficient functions. These expressions 
serve as a check on the exact calculation in [Q which is only available in a 
large computer program involving numerical integrations over several vari- 
ables. Furthermore it enables us to see at which Q^-values the asymptotic 
heavy quark coefficient functions coincide numerically with the exact ones, 
which gives an indication when the charm quark can be treated as a massless 
quark. 

The content of this paper can be summarized as follows. In section 2 we 
introduce our notations and give an outline how the heavy quark coefficient 
functions can be determined in the asymptotic limit ^ m^. In section 3 
we present the calculation of the full two-loop operator matrix element needed 
for the computation of the asymptotic form of the heavy quark coefficient 
function. The latter will be calculated in section 4. Finally in section 5 we 
will show at which Q^-values the asymptotic form coincides with the exact 
one given in 0]. In Appendix A an exact analytic expression for the heavy 
quark coefficient function, which is valid for any and m^, is presented in 
the case of the Compton subprocess. An important trick how to compute an 
operator matrix element with five different propagators is given in Appendix 
B. The long formulae obtained for the full operator matrix elements and the 
asymptotic heavy quark coefficient functions are presented in Appendices C 
and D respectively. 
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2 Heavy flavour coefficient functions 



In the section we will show the connection between the heavy flavour coeffi- 
cient functions computed in the asymptotic limit ^ and the operator 
matrix elements (OME's). Here denotes the mass squared of the virtual 
photon with momentum q [q^ = —Q^ < 0) and m stands for the heavy 
flavour mass. The variable x is deflned by x — Q^/2p ■ q (Bjorken scaling 
variable) where p stands for the momentum of the proton. The OME's arise 
when the local operators, which show up in the operator product expansion 
(OPE) of two electromagnetic currents, are sandwiched between the proton 
states indicated by |P >. In the hmit that x is fixed and ^ M^, where 
M denotes the mass of the proton, this OPE dominates the integrand of the 
hadronic structure tensor W/^u, which is defined by 

{p,q) = j^J d'ze^^-^ <P\[J,{z),MmP> ■ (2-1) 

This structure tensor arises when one computes the cross section for deep 
inelastic electroproduction of heavy flavours 

e-(£i) +P(p) ^ e-(4) + g(pi)( Q{P2) )+'X', (2.2) 

provided the above reaction is completely inclusive with respect to the hadronic 
state 'X' as well as the heavy flavours Q{ Q )■ When the virtuality { q = 
ii — ^2 ) of the exchanged photon is not too large ( <^ M| ) the reaction 
in (2.2) is dominated by the one-photon exchange mechanism and we can 
neglect any weak effects caused by the exchange of the Z-boson. In this case 
Wni,{p, q) in (2.1) can be written as 

W^,u{p, q) = -^{a^u - '^^)Fl{x, Q^) + {p^,Pu - ^{p^qu + Puq^) 

^ p.q 

The above formula follows from Lorentz covariance of Wfj_i,{p, q), parity invari- 
ance and the conservation of the electromagnetic current J^. Besides the lon- 
gitudinal structure function Fl{x,Q'^) and the structure function F2{x,Q'^) 
one can also deflne the transverse structure function Fi{x,Q'^) which, how- 
ever, depends on the two previous ones and is given by 

F,{x,Q') = ^[F2{x,Q') - Fl{x,Q')] . (2.4) 
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The above structure functions show up in the deep inelastic cross section of 
the process in (2.2) provided we integrate over the whole final state 

^ - l^^'IO + (1 - yrm^.Q') - y'F^M')] . (2.5) 

where S denotes the square of the cm. energy of the electron proton system 
and the variables x and y are defined ( see above (2.1) ) as 

x = ^ (0<x<l) , (0<y<l), (2.6) 

2p ■ q p ■ ti 

with 

- g2 = Q2 ^ xyS . (2.7) 

In the QCD improved parton model the heavy flavour contribution to the 
hadronic structure functions can be expressed as integrals over the partonic 
scaling variable z = Q"^/ (s + Q'^), where s is the square of the photon-parton 
centre-of-niass energy (s > 4m^). This yields the following result 

F,{x, g^ m') = X -[-T. el{^{^, ^^')lUz, Q\ m\ /.^) 

/.^)l,,(., Q\ m\ /.^)} + A(^, ii^)L^l{z, Q\ m\ /.^)' 

+a; e\ f"'"' — {e(-, ix')E,^^{z, Q^ m\ ix^) 

Jx Z ^ ^Z ' 

+G(^,^2)i/,,(^,g^m^/.2)}, (2.8) 

where i = 2,L and the upper boundary of the integration is given by Zmax = 
Q^/{Am'^ + Q^). The function G{z,iJ,^) stands for the gluon density. The 
singlet combination of the quark densities is defined by 

n f 

^z,^l')^Y.iM^^^^")+7^{z,^l')), (2.9) 

i=l 

where fi and /j stand for the hght quark and anti-quark densities of species 
i respectively. The non-singlet combination of the quark densities is given by 

A{z, /.^) = E - - E 4) {M^, i^') + f^')) ■ (2.10) 

i=i k=i 
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In the above expressions the charges of the hght quark and the heavy quark 
are denoted by and en respectively. Furthermore, Uf stands for the number 
of hght quarks and /i denotes the mass factorization scale, which we choose 
to be equal to the renormalization scale. The latter shows up in the running 
coupling constant defined by as{fi'^). Like the parton densities the heavy 
quark coefficient functions Lij {i = 2, L; j = q, g) can also be divided into 
singlet and non-singlet parts which are indicated by the superscripts S and 
NSineq. (2.8). 

The distinction between the heavy quark coefficient functions Lij and 
Hi j can be traced back to the different photon-parton production procesess 
from which they originate. The functions Lij are attributed to the reactions 
where the virtual photon couples to the light quark, whereas Hij originates 
from the reactions where the virtual photon couples to the heavy quark. 
This explains why there is only a singlet part for Hi j and that Lij, Hi j 
in eq. (2.8) are multiplied by e^, respectively. In P the heavy quark 
coefficient functions have been calculated up to next-to- leading order (NLO). 
In the Born approximation (first order of as or LO) one has the photon-gluon 
fusion process 

i*{q)+9{h)^QiPi) + Q{P2), (2.11) 

which leads to the lowest order contribution to Hi g denoted by H^^g . The 
next order is obtained by including the virtual gluon corrections to process 
(2.11) and the gluon bremsstrahlung process 

1*{q)+9{ki) ^ 9ih) + Q{pi) + Q{p2) , (2.12) 

(2) 

both of which contribute to the second order term in Hi^g denoted by Hlg . In 
addition to the above reaction we also have the subprocess where the gluon 
in (2.12) is replaced by a light (anti-) quark, i. e. 

7*{q) + qiq){h) qiq^h) + QiPi) + Q{P2) ■ (2.13) 

This process has however two different production mechanisms. The first one 
is given by the Bethe-Heitler process (see figs. 5a,b in f^) and the second 
one can be attributed to the Compton reaction (see figs. 5c, d in [Q]). In 
the case of the Bethe-Heitler process the virtual photon couples to the heavy 
quark and therefore this reaction contributes to Hig. This second order 
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(2) 

contribution will be denoted by H-^J . In the Compton reaction the virtual 
photon couples to the light (anti-) quark and its contribution to Lf^^ will be 

denoted by L^^'^^''. Since can be written as Lf^^ = Lf^ + L^^, L^^'^^^ 
also contributes to the singlet part of the coefficient function. In general the 
heavy quark coefficient functions are expanded in as as follows 



i^,,(^,Q^m^/i2) = (^) v;(^,Q^m^/i2) , (2.14) 



i^,,(.,Q^m^/.^) = E {^yHt^iz,Q^m\f,') , (2.15) 

k=2 ^'^ 



00 

^/z,Q^m^^,') = Y.{^) Llf\z,Q\m^^,^), (2.16) 

k=2 

where r = NS,S. Finally we want to make the remark that there are no in- 
terference terms between the Bethe-Heitler and Compton reactions in (2.13) 
if one integrates over all final state momenta. 

The complexity of the second order heavy quark coefficient functions pro- 

NS (2,) 

hibits publishing them in an analytic form, except for L^^ , which is given 
in Appendix A, so that they are only available in large computer programs 
involving two-dimensional integrations. To shorten the long running time 
for these programs we have previously tabulated the coefficient functions in 
the form of a two dimensional array in the variables rj and ^ in a different 
computer program [jl2[. These variables are defined by ( using s = {q + ki^, 
z = Q'^/{s + Q^), and r] = {s - Am^)/{Am^)) 

V=^-^^-l , e = ^. (2.17) 

This new program has shortened the computation of the charm structure 
functions Fi{x , Q"^ , 771"^) considerably as one only requires one integral over 
the variable z in (2.8), therefore making our results for the NLO corrections 
more amenable for phenomenological applications. However, when 3> 
it is possible to get complete analytic forms for the heavy quark coefficient 
functions which are similar to the ones presented for the light quark and 
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gluon coefficient functions given in |T^. To get the analytic form in the 
above asymptotic regime one can follow two approaches. The first one is 
to go back to the original calculation of the exact coefficient functions in 
[0] and repeat the computation of the Feynman graphs and the phase space 
integrals in the limit ^ iv? . An example of such a calculation can be 



found in [|T3|, where all photonic corrections to the initial state of the process 
e" + e"*" ^ /i~ + /i"*" in the limit S 3> ml were computed. However, this 
procedure is still quite complicated because one cannot neglect the fermion 
masses at too premature a stage which results in rather messy calculations. 
Fortunately, as one can find in [|1^], there exists an alternative method which 
we will use for the heavy quark coefficient functions. 

In the limit ^ m? the heavy quark coefficient functions behave loga- 
rithmically as 

= E4''(-,4)in'g , (2.18) 

^ -j^ I I L lib 

with a similar expression for L-^-^ As has been already mentioned in the 
introduction these large logarithms \vl{Q'^ /w?) dominate the radiative cor- 
rections. This is in particular the case for charm production in the large 
region which is accesible to HERA experiments. The above large logarithms 
also arise when is kept fixed and m? so that they originate from 
collinear singularities. These coUinear divergences can be removed via mass 
factorization. The latter proceeds in the following way. In the non-singlet 
case we have 

G-(^. „,) + ^) = r-(^) » n, + 1) , (2.19) 

For the singlet case the mass factorization becomes 

T.n^{^)^Cl&nj + l), (2.20) 

with k,l = q, g. Notice that in the above expression we have suppressed the 
z-dependence for simplicity. The convolution symbol is defined by 

(/ ® g) (z) = £ dzi £ dz2 5{z - z^Z2)f{z^)g{z2) . (2.21) 
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The quantities Tki and Ci^k {k, I = <l-,g) which appear in the above equation 
stand for the transition functions and the hght parton coefficient functions 
respectively. Notice that all mass dependence is transfered to the transition 
function F^/. The removal of the logarithmic terms ln'(/i^/m^) from the heavy 
flavour coefficient functions leads to an enhancement by one of the number of 
light flavours Uf in the light parton coefficient functions Ci^k - see eqs.(2.19), 
(2.20). The latter have been calculated up to second order in in Since 
these coefficient functions do not depend on the way one has regularized the 
collinear divergences in the parton cross sections (for a discussion see fl^) 
one can also use them in equations (2.19), (2.20). Now if one also knows the 
transition functions Tki one can reconstruct the asymptotic behavior of the 
heavy flavour coefficient functions Li^k and Hi^^. This is possible because the 
transition functions T^i are identical to the operator matrix elements (OME) 
denoted by A^i. The latter appear in the operator product expansion (OPE) 
of the commutator of the two electromagnetic currents in (2.1) near the 
lightcone. Suppressing the Lorentz index of the electromagnetic current the 
expansion can be written as 

hm [J(z), J(0)] =Y.Y.Ct\^'W. ■ --^.^Or'-iO) , (2.22) 

k m 

where the distributions C^\z'^) {k = q,g) are the Fourier transforms of 
the light parton coefficient functions C^\Q'^ / (i'^) defined in (2.19), (2.20). 
Notice that the latter has been Mellin transformed according to 

In the above expressions the index i [i = 2, L) in C^^^ has been omitted since 
we suppressed the Lorentz index of the current J^{z). Inclusion of the latter 
index implies that one has two independent structure functions F2 and Fl. 
The quantity m refers to the spin of the local operators 0^1^'"^"" which appear 
in the OPE (2.22). The latter are given by the non-singlet quark operator 
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A. 



0^j,-'""(x) = -i'^-^SYl>{x)-i^'^D^'^ ■ ■ ■ D^'"' ^i){x)\ + trace terms, (2.24) 



and the singlet operators 



i;{x)Y'^D''^ ■ ■ ■ D^'^ipix) + trace terms , (2.25) 
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= li'^-^s\F^''''{x)D^' ■ . •D''— + trace terms . 

(2.26) 



Here S denotes the symmetrization of the operators in their Lorentz indices 
ijLi and the trace terms are needed to make them traceless. The in (2.24) 
represent the generators of the flavour algebra, whereas the index a in (2.26) 
is the colour index. The objects ip{x), represent the quark field and the 
gluon field tensor respectively and denotes the covariant derivative. 

The OPE expansion in (2.22) can be applied in the limit 3> and 
fixed X so that the integrand in (2.1) gets its dominant contribution from 
the lightcone behavior of the current-current commutator. Inserting (2.22) 
in (2.1) and replacing the hadron state |P > by a light quark or gluon one 
gets the relations 



, (2.28) 



where i = 2, L; k,l — q, g. Since the operators in (2.22) are already rcnormal- 
izcd, the above expressions can be only collinearly divergent. The collinear 
divergences due to the presence of light partons are regularized by the method 
of n-dimensional regularization. They are indicated by the pole terms (1/e)* 
(e = n — 4). However the collinear divergences due to the heavy quarks 
are regularized by the mass m which shows up in the form of logarithms of 
the type \n\Q'^/m'^), W {fi'^ / m^) . The objects jFj^^ denote the light parton 
structure functions. The OME's A^i are defined by 

=<W> , (2.29) 

where |Z > is a light quark or a gluon state and all quantities depending 
on the pole terms are indicated by a tilde. If the latter are removed via 
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mass factorization one obtains expressions (2.19), (2.20) which imphes the 
following identification 



r?^(^) = <,^(^). (2-30) 



For the computation of the asymptotic behaviour of the heavy quark coef- 
ficient functions corresponding to the processes (2.11)-(2.13) one needs the 
following quantities. For processes (2.11) and (2.12) we have to calculate the 
one-loop OME Ag^g and the two-loop OME Agg respectively. They are repre- 
sented by the Fcynman graphs in fig.l and fig. 2. The Bethc-Hcitlcr process 

PS f2l 

given by reaction (2.13) corresponds to the two-loop OME Ag^ in fig. 3 
whereas the Compton process (2.13) is related to A^^'^'^^ with the two-loop 
OME's in fig. 4. The calculation of the graphs in figs. 1-4 and the derivation 
of the asymptotic form of the heavy quark coefficient functions will be the 
aim of the next two sections. 
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3 Calculation of the two-loop operator ma- 
trix elements 



Before presenting the results of our calculation of the OME's we will first 
derive the general structure of the OME's discussed in the last section. If 
we insert the OPE (2.22) into the structure tensor PV^^ in (2.1) the OME's 
which are derived from the Feynman graphs in figs. 1-4 will be computed in 
the forward direction. The latter means that the momentum leaving the 
operator vertex equals zero. Further if one puts the momentum, indicated 
by p, of the external light quark and gluon off-shell (p^ < 0) only ultraviolet 
(UV) singularities appear in the OME's. Using this off-mass-shell assignment 
one can express the renormalized as well as the unrenormalized OME's into 



the renormalization group coefficients as is done in |T^. However for our 
computations we have to put the external momentum on-shell {p^ = 0) so 
that the OME's turn into genuine S-matrix elements. This mass assignment 
implies that in addition to UV divergences one also encounters collinear (C) 
divergences which originate from the coupling of the external on-shell mass- 
less quanta to internal massless quanta. In the computation of the Feynman 
graphs both types of divergences will be regularized using the technique of 
n-dimensional regularization. However since both singularities in the OME's 
will manifest themselves in the form of pole terms of the type e"'^ (e = n — 4) 
it is very hard to trace back their origins. Nevertheless one can express the 
OME's into the renormalization group coefficients in a similar way as has 



been derived for the off-shell case in Where possible we will make a 

distinction between UV-pole terms and C-pole terms, which are indicated by 
e^Y and e^^ respectively , and identify them euv = when it is appropriate. 

In the subsequent part of this section we will construct the OME's, cor- 
responding to the graphs in figs. 1-4, in such a way that the coefficients of 
the pole terms are given by the renormalization group. These coefficients are 
products of the terms appearing in the beta-function and the AP-splitting 
functions (anomalous dimensions) [|1^. The purpose of this presentation is 
threefold. First we need these coefficients for the construction of the heavy 
quark coefficient functions in the next section. Second since the renormaliza- 
tion group coefficients are known in the literature we can predict the residues 
of the pole terms so that these expressions serve as a check on our calcula- 
tions. Third it is much easier to show the renormalization and mass factor- 
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ization for the algebraic expressions, which are short, than for the analytic 
formulae in our calculations because the latter are rather long. 

The OME's can be expanded in a perturbation series as follows 

^« = £©X"- (3.1) 

fc=0 ^" 

In the following discussion we distinguish three different types of OME's. 
First we have the unrenormalized ones indicated by Aij. They contain UV- 
as well as C-singularities. Second after renormahzation the UV-divergences 
are removed and we are left by the OME's defined by Aij which still contain 
C-divergences. The latter have to be removed via mass factorization so that 
the Aij turn into the finite OME's indicated by Aij. Notice that the expansion 
in (3.1) holds for all three different types of OME's. 

The renormahzation of Aij proceeds in three steps. First we will perform 
mass renormahzation for which we choose the on-mass-shcll scheme. This 
implies that the bare mass m, which occurs in Aij has to be replaced by 

rh = mil + ^6m) , 6m = CfSJ^Y^^ \— - a} , (3.2) 



so that the mass renormalized Aij reads up to order 

4 = 4? + (|)4' + (S)Vl? + *'»5l4"}- (3.3) 

Notice that the zeroth order term A^'^^ = 6ij is mass independent. In 
the above perturbation series the quantities //^ and are artefacts of n- 
dimensional regularization. The mass parameter n originates from the di- 
mensionality of the gauge coupling constant g (a^ = g'^/An) in n-dimensions 
and should not be confused with the renormahzation and mass-factorization 
scale. However if one only subtracts the pole-terms like in the MS-scheme 
the mass parameter /i turns into the afore-mentioned scales. The spherical 
factor is defined as 

5, = exp{i(7E-ln47r)}, (3.4) 
and the colour factor is 

- 1 
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in SU(N). Further ctg denotes the bare couphng constant which will be renor- 
malized as follows 



s 



a,=aV)(l + ^^'5a.), (3.6) 
S.{^ + ^(4)"''(l + 5<(2)^w)}- 



Here a summation over the heavy quarks H{H = c,b,t) is understood and 
C(2) = TT^/ 6 is the Riemann zeta- function. Further Pq is the lowest order term 
in the series expansion of the beta-function, which up to two-loop order, is 
given by 

m--^.£,-P.jSf. (3.8) 



where 



and 
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M 20 
Pi = y C-l - iCpTfUf - -C/TfUf . 

Here Ca and Tf denote colour factors of SU(N) 

Ca = N , Tf = \, (3.9) 

and Uf denotes the number of light flavours which enter via the fermion- 
loop contributions to the gluon self-energy. However besides the light quarks 
also the heavy quarks with mass itlh contribute to the renormalized coupling 
constant a^. This contribution is indicated in (3.7) by I3q_h = —4:Tf/3 and 
niH > TTi. The coupling constant renormalization in (3.6), (3.7) is determined 
in the MS-scheme as far as the light flavours and the gluon are concerned. In 
addition we make the choice that the heavy quarks decouple in the running 
strong coupling constant asi^"^) for /i^ < mjj and the renormalized OME's. 
The factor l+e^C{2)/8 in (3.7) arises from the requirement that H//(0, m^) = 
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0, where Unip'^j'm'jj) is the contribution to the gluon self-energy due to the 
heavy quark loops indicated by H. After coupling constant renormalization 
Aij takes the form up to 0{a'^) 

The remaining UV-divergences are removed by operator renormalization 
which is achieved by 

Ai(euv,ec, = ^ifc(euv, a^) ® ^fci(ec, "s) , (3-11) 

where Zij = q, g) are the operator renormalization constants correspond- 
ing to the operators in (2.24)-(2.26). Notice that for the non-singlet operator 
in (2.24) Z^^ is a real number, whereas for the singlet operators in (2.25), 
(2.26) Zij becomes a matrix. The operator matrix elements can be expanded 
in as analogous to A^j in (3.1) as follows 

oo , 

% = E(^)4'. (3-12) 

so that the renormalized OME's Aij in (3.11) read up to 0(af), (using Zij^ = 

+fo,iW + Z-,«'>ig' + Z-/-''>} , (3.13) 

where Zj^^''"'^^ are the expansion coefficients of the inverse matrix Z^^. Choos- 
ing the MS-scheme one obtains the following expression (see [|15]) up to 0(af) 

Z^j{evy,as) = S^j + (^)Se[ —Pip] 

^ 477 ^ euv 




(3.14) 



Here P}'j'^ and Pj^j^ stand for the first and second order AP-splitting func- 



tions |16] respectively which have been calculated in the literature (see [It 

|T8|, ||20|). Notice that in our notation Aij and depend on the 
partonic Bjorken scaling variable z in (2.8) which for convenience has been 
suppressed in our formulae. Hence instead of multiplications we have to deal 
with convolutions denoted by the symbol ®, which is defined in (2.21). Fur- 
thermore the AP-splitting functions are related to the anomalous dimensions 
7!™'* of the composite operators in (2.24)-(2.26) via the Mellin transform 



(m),(fe) 
Jij 



dz z^'-^P.fiz) 



(3.15) 



Finally we have to remove the C-divergences appearing in (3.11). This 
is achieved by performing mass factorization which proceeds as 



a. 



(3.16) 



where Aij denote the finite OME's which do not have UV or C divergences 
anymore. The quantities Tij stand for the transition functions which have 
the same properties as mentioned above (3.12). Like Aij and Zij the Fjj 
can be expanded as a power series in 



k=0 



(3.17) 



so that the mass factorized OME's Aij read up to O(a^) (using F^-^-* = 6ij] 
Ai 



Y 

47r/ 



+{i 



d 



i(i) 



(1) 



'-1,(1) 4(1) 



^-1,(2) 



(1) 
ik 



dm 

-1,(1) 



ik 



kj 



f:" 



-1,(2) 



(3.18) 



where F,/'*-'^'* are the expansion coefficients of the inverse matrix F ^. Choos- 



ing the MS-scheme one obtains the following expression (see [|TT]) up to 0(q;^) 



r.,(ec,a.)=5,, + 0Sjlpf 
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+{^fsnir{lp!^ » PUf + fti^f } + ■ (3.19) 

If all quarks would be massless we would have the identity Tij = . How- 
ever since the heavy quark Q is massive it does not contribute to those 
splitting functions P^^ which appear in the transition functions. Hence in 
this case i and j only represent the light quarks and the gluon. This is in 
contrast to the operator renormalization constant Zij where i and j can also 
stand for the heavy quark. The same assertion holds for the heavy flavour 
contributions to the beta-function given by /3o,i? — Q^) which shows up 
in (3.14) but not in (3.19). The reason is that the mass m of the heavy quark 
acts as a regulator for the C-divergences but not for the UV-singularities. 

Using the master formula in (3.18) one can now construct the general form 
of the finite OME Aij as well as the unrenormalized OME Aij expressed in 
their renormalization group coefficients /9o, f3o,H and P^j^ {k = 0,1; i,j = 
q, g). For convenience we will take for A^j the representation (3.3), where the 
mass renormalization has been already carried out. 

Starting with the one- loop OME A^q], which receives contributions from 
the graphs in fig.l, we obtain 



4] - - -^Pi? + 41 + ^waS} . (3.20) 



A*^ ' ' euv 
The renormalization group coefficients are given by 

P(;) = STf[z' + (1 - zf] , 

«S = -^C(2)Pi?. (3.21) 

The expression for the two- loop contribution to Aqg (see fig. 2) can be written 
as 
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+ 



(1) 

Qg ■ 



(3.22) 



Notice that in the above expression the pole terms e'^ stand for the UV as 
well as C- divergences so that we have put euv = ec- Further we infer from 
the literature that 



Pr=4C.[2(^ 

= ^^^ii-r—)^ + --2 + z-z']+ 2/3o5(l - z) , 
p« = 8CFTf[2{l - 2z + 2z'^){\n\l - z) - 21n;2ln(l - z) - 2C(2)} 
+ (1 - 2;z + 4^2) In^ z + 8;z(l - z) ln(l - z) + {?, - Az + 8^^) Inz 



+ 14 - 29^ + 2Qz' 



+8CATf [(1 + 2^ + 2z'^){\n'^ z-A Inz ln(l + - 4Li2(-2) - 2C(2)} 
+2(1 -2z + 2z^) [C(2) - ln2(l - z)] - (3 + 6^ + 2z^) In^ z 



88 

-8z(l - z) ln(l - 2) + (2 + 16z + — z^) In 2 



+-f— - 18 + 225^-21^ 
9^2; 



(3.23) 



and Oq] (3.22) has to be computed in this paper. In the above and subsequent 
expressions the functions Li„(2;) denote the polylogarithms which can be 
found in PT|. The finite OME's follow from (3.18) where we have 



aW _ ^(1) , 7-1,(1) 



(3.24) 
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and 



M _ 2W I i„ 4(1) I z-^M a 2W + z-i.(i) a 



,(1) 



+z,-M^i + (iS + z-M'))»r-' 

+{ - i/^,') - Aag - i/^J) » ag + ip<,°) » ag} In 



(3.25) 



In eqs.(3.24), (3.25) wc have used the heavy flavour contributions (m|^ > Q^) 
to the one loop OME in fig.lc which is given by 



2 1 

Po,H - 7euv/5o,//C(2) 

euv 4 



(3.26) 



and the constants 



P(0) 
9ff 



Cuv "-2 



53 



2euv 
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;7-i,(i) _ c rJ_p(o) 
«^ ^Uuv '''' 



55 



55 



(3.27) 



Notice that due to our scheme for Us in (3.7) /3o,// and m// (m^ > Q^) have 
completely disappeared from (3.25). 
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^PS (2) 

The renormalization and the mass factorization of the OME's Aq^ in 
fig. 3 proceeds in the same way as done for Agg. The superscript PS stands 
for 'pure singlet' and it originates from the definition that the singlet OME 
Agg can be decomposed into 



A 



PS 



^qq +^qq 



(3.28) 



where A^^ stands for the non-singlet OME. The unrenormalized OME reads 



PS,(2) 



(3.29) 



Like in the case of Aq^ we did not make any distinction between UV and 
C-singular pole terms e^'^ (euv = ec)- The renormalization group coefficients 
are given by (see also (3.21)) 



pPS,(i) ^ g^^^jy [ - 2(1 + z) In' ^ + (2 + 10^ + ^z^) In 



40 , _ 112 2 



(3.30) 



and agg' will be calculated in this paper. The finite OME Aq^^^' can be 
again derived from (3.18) which yields 

-1,(2) 



4PS,{2) 



1PS,(2) 



+{-k?''"+k>^':'}"> 



a9? + z-^'W^ (s?) r-^'(^) 



m 



1^' 



m 



1^' 



PS,(2) 



(3.31) 
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Here we have used the renormahzation group constants (see also (3. 21), (3. 30)) 



-1,(2) 



L/ip(o) 

i \o 19 



p(0)\ I ^ pPS,(l) 

J 2euv 
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J_p(0) 

ec 



(3.32) 



The last OME which we have to deal with is represented by the heavy 
quark loop contribution to A^^'^'^^ in fig.4. The unrenormalized expression 
reads 



JNS,(2) 

^qq,Q 



2 

m \e 



'■{-/»o,«<'} + — {-;pS"} + <': 



-uv 



1 ^NS,(1)-, , NS,(2) 

2 - qq,Q i ^qq,Q 



(3.33) 



Contrary to and ^q^'^^^ the above OME only contains UV-divergences 

since the heavy quark Q prevents -^q^q^ to be C-singular provided we choose 
the coupling constant renormahzation scheme in (3.7). The renormahzation 
group coefficients in (3.33) are given by 



l3o.. 



Q 



,NS,(1) 



CfTj 



160 



9 VI _ ^ 



176 16 161 + ^2 

H z In z 

9 9 3 1-z 



+5(l-.)(-^-fc(2) 



(3.34) 



Since A^^^'^^^ does not contain C-divergences we only need operator renor- 
mahzation to render it finite. From (3.18) we infer 



.NS,(2) _ iNS,(2) , 

^qq,Q - ^qq ^ \^qq,Q 



(3.35) 
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with 



-uv ^uv M 

. (3.36) 



Notice that expressions (3.33) and (3.35) also contribute via (3.28) to the 
singlet part of the OME. 

Since we can infer the coefficients of the double and single pole terms 
of the unrenormalized OME Aij from the two-loop corrected AP-splitting 
functions and the beta-function, the above expressions serve as a check on our 
calculations. The non-pole terms defined by a\f , which cannot be predicted, 
have to be calculated in this paper. They are needed to compute the heavy 
quark coefficient functions (2.18) up to the non- logarithmic term, which will 
be done in the next section. 

Before finishing this section we will give an outline of our calculation of the 
OME's depicted in figs. 1-4. We have computed the OME using the standard 
QCD Feynman rules and the operator vertices corresponding to (2.24), (2.25) 
which can be found in the literature (see [1^). Since in our case the latter 
are S-matrix elements we have to consider the connected Green's functions 
where the external quark and gluon propagators are amputated. Notice that 
one has to include the external self energies. The connected Green's function 
needed for Agg (figs. 1,2) is given by 

< 0|T(A;:(a;)Of ■••'^-(0)4(2/)|0 (3.37) 
and for A^^ (fig-3) and A^^ (fig-4) we have 

< 0|T(V'.(x)O^^'^-^'-(0)^,(i/)|0 (3.38) 

with r = PS and r = NS respectively. Further a, b and i,j are the colour in- 
dices of the gluon field and the quark field ip respectively. Since the above 
operators O'^^""'^'"(0) are traceless symmetric tensors under the Lorentz group 
the computation of the connected Green's functions reveals the presence of 
many trace terms which are not essential for the determination of the anoma- 
lous dimensions or splitting functions. Therefore we will project them out by 
multiplying the operators by an external source J^^....^^ = A^j...A^,,^ with 
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= 0. Performing the Fourier transform into momentum space and sand- 
wiching the connected Green's functions by the external gluon polarizations 
and quark spinors one obtains 

e''{p)Gt.e''{p) , (3.39) 

and 

u{p,s)Gl^Xru{p,s), (3.40) 

where \r denote the generators of the flavour group SU(n/). The tensors 
G"''„„ and Gi^ have the form 

Gt. = ^ir^ {e, ^, a,)6^'{A • p)-( - g,^ + ^^^^^^) , (3.41) 

and 

^A,,{e,%,aM\A-pr-'4. (3.42) 
By projecting these tensors out one obtains finally the OME's 

- ^,:^2^-9.^)nA-p)-Gt,., (3.43) 



^'"^) = . pr'^Tri^G^^) . (3.44) 

Notice that in (3.43) the summation over the dummy indices /i and v includes 
unphysical (non-transverse) gluon polarizations which have to be compen- 
sated by adding graphs containing the external ghost lines represented in 
fig.2s and fig.2t. Instead of —g^^'^/{n — 2) one can also use the physical po- 
larization sum [—Qnv + {^uPv + Ai,p^)/A • p\/{n — 2) and omit the graphs 
with the external ghost lines. However in this case the individual Feynman 
graphs lead to integrals with higher powers of the momenta in the numera- 
tor so that they become more difficult to compute. Moreover the number of 
independent integrals is artificially increased. The advantage of constructing 
the Green's function in the way shown in (3.43), (3.44) is that one does not 
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have to resort to complicated tensorial reduction programs as had to be used 
for example in [ 1^ . Therefore one has only Lorentz scalars in the numerators 
of the Feynman integrals which can be partially cancelled by terms in the 
denominators. The traces of the fermion loops in figs. 1,2 and the contraction 
with the metric tensor in (3.43) have been performed by using the algebraic 
manipulation program FORM [^. We did the same for the graphs in figs. 3,4 
where we had to compute the trace in (3.44). The computation of the scalar 
integrals is straightforward as long as the number of propagators does not 
exceed four. In the case of five different propagators the calculation becomes 
more cumbersome but here one can use the trick of integration by parts 
2^ . Examples are given in Appendix B. The results for the unrenormalized 



OME's are too long to put them in this section and we will defer them to 
Appendix C. 

Apart from the check on the pole terms mentioned above we can also check 
the finite term of the Abelian part of Aqg. Removing the overall constant 
(1 + (— 1)'")/2 the Mellin transform reads up to order (see Appendix C) 



7 



i(2),(m) ^ J(2),(m)- 

(3.45) 



in the above equation so that we get rid of the term proportional to in 
(C.l). If we now take the first moment m = 1, Aq^'^^^ ^ (m — 1)^^ but 



Notice that coupling constant renormalization has already been carried out 

f 

A-Q^g^^^ and Aq^g^P are finite. One can now easily show that there exists a 
relation between the Abelian terms and the heavy fermion loop contribution 
to the gluon self energy which we will denote by 

^v(p^m^) = S,(^)^,^«(p^m^) + S,^(^)UT,^^,^^(p^m2) . (3.46) 

The relation is given by 

n?)(0,m^) = liS«, (3.47) 



U^\0,m^) = lA^lf. (3.48) 
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flS^^^(0,m2) and flS,^^(0, m^) can be inferred from eqs. (2.11) and (5.1) in 
In the latter reference the heavy fermion loop contribution to the self energy 
of the Z-boson and the photon were calculated. The self-energy contribution 
to the gluon and the photon are related as follows 

fI^,^)(p^m2)=4l^^(p2), (3.49) 



2 2\ 

■ ,m ) 



(3.50) 



where Ilj'{p'^) in eqs. (3. 49), (3.50) are equal to the expressions quoted in 
(2.11) and (5.1) of respectively. The factors of 4 originate from the 
convention that in one expands the quantities in a^/vr instead of as/ (47r) 
as is done in this paper. The factors of 4/3 and 3 refer to the colour factors 
Cf and Ca in SU(3). From |^ we infer that 





r 8n 




~ 3e 



(3.51) 



+ 15 



(3.52) 



where n = 4 + e. (Notice that in |2J] n = A — 2e). Calculation of the first 
moment of (3.45) (see also (C.1),(C.2)) reveals that the relations in (3.47), 
(3.48) are satisfied. 
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4 Heavy quark coefficient functions 



In this section we will compute the heavy quark coefficient functions Hi^i and 
Li^i defined above (2.14) up to order in the asymptotic limit ^ m^. For 
this computation we will use the mass factorization theorems as represented 
by eqs. (2.19) and (2.20) where the transition functions T kii^Ji^ / m?) stand 
for the finite OME's Aki{ix^ /m?) computed in section 3. 

Let us start with the heavy quark coefficient functions which originate 
from the subprocesses where the virtual photon couples to the heavy quark. 
Here the mass factorization theorem implies 

where Hi^i [i = 2,L]l = q, g) denote the heavy quark coefficient functions in 
the limit ^ and Am {k, I = q,g) are the finite OME's computed in 
the last section. The coefficient functions Ci^i{Q'^ / fi^) have been calculated 
ri| . They are obtained from the massless parton structure functions as 



m 

defined in eqs. (2.27), (2.28), after having performed mass factorization in 
the MS-scheme. In the last section also the OME's A^i have been calculated 
in the MS-scheme. In the product on the right-hand-side of (4.1) the scheme 
dependence is only partially cancelled which is revealed by the fact that Hi i 
is still scheme dependent. This dependence, which is indicated by ix^/w? 
in Hi I , originates from the coupling of a light parton (gluon or quark) to 
an internal light parton characteristic of the production mechanisms of the 
processes in (2.12) and (2.13). 

In the case of the lowest order photon-gluon fusion process we have (see 
(2.11) and (2.14)) 

<{-.A)=cU^). (4.2) 



A^M( 



(4.3) 



where 

are the order gluonic coefficient functions which e.g. can 



be found in |Tl| . Like the OME's they can be expressed in the renormalization 
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group coefficients as 



cS(f)=cS. (4.4) 



From (3.24) and (4.2)-(4.5) one obtains the asymptotic forms for the order 
tts heavy quark coefficient functions 

<(S.^)=c«, (4.6) 



< ^) = \pS^ In % + + eg . (4.7) 



In order al the coefficients of the photon-gluon fusion process (2.12) become 



<(S.£;)-4^(£;)«cs(f)+cag). (4.8) 



^i3(S. ^) - + « ) . c-S(f ) , (4.9) 

where C\]^, (« = 2, L) become (see JUl) 

cS(f)=d'.',. (4.10) 

+4:^ (4.11) 

(4.12) 



29 



+42 . (4.13) 

Prom (3.24), (3.25) and (4.10)-(4.13) one infers the asymptotic form 

(4.14) 



„2 ■ L4-I9 --59 2""" '» ' ^ ^ (1= 
+{ - /3o(4:J + '■S) + ^^-r ® (4'] + «S)} In ^ 

(4.15) 

Notice that the above formulae are still dependent on the mass factorization 
and renormalization scale yU^. The same scale dependence was also found 
for the exact expressions for the heavy quark coefficient functions where 
and VP? can be arbitrarily chosen. The /x^ dependence can be attributed to 
coupling constant renormalization represented by the lowest order coefficient 
/3o in the beta-function and the lowest order splitting function P'^^ standing 
for the transition g ^ g + g. 

The computation of the asymptotic expression of the heavy quark coeffi- 
cient function corresponding to the Bethe-Heitler process (2.13) proceeds in 
the same way. From (4.1) we derive 
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2,9 V77i2 ' jyj2 



(4.17) 



The coefficient functions Cf^q'^'^^ {i 
read 



2, L) are computed in |TT|] and they 

(f ) = {^nT ® -S} ^ + ^J'" . (4.18) 



(4.19) 



, PS,{2) 



Using (3.31) and (4.18), (4.19) the heavy quark coefficient functions become 



- « 1" S + {in'?' « c« } .„ ^ 



+C 



PS,(2) 
L,q 



(4.20) 



+cr,'^'+«S''''-nT®<'S- (4.21) 

Notice that the above /i^-dependence can be only attributed to the mass 
factorization scheme which enters via the transition q ^ g + q represented 
by the sphtting function Pgg^- The mass factorization for the heavy quark 
coefficient functions where the photon couples to the light quark proceeds in 
a similar way as given in eq.(4.1) i.e., 



2 j ~ "^^'l 2 



(4.22) 
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In the case of the Compton process (2.13) where Li^q is of order one can 
derive 

^L,, - [-^) , (4.23) 



-NS,i 



where Q^^j (^ = 2, L) denotes the quark loop contribution to the non- 
singlet coefficient function for one specific massless flavour / which is given 
by (see H) 



^Ns,(2)/QS_r In p(o)\i^2 , rlpNs,(i) ^ 



+cSf . (4.26) 

Notice that the above coefficient functions are represented in the MS-scheme 
for the coupling constant renormalization (see (3.6)). However A^^'^"^^ in 
(3.35) has been determined in another scheme where the heavy flavour Q 
decouples in the strong coupling constant (see (3.7)). Choosing the latter 
scheme and putting f = Q we obtain instead of (4.25), (4.26) the expressions 

with Cl^ deflned in (4.10), (4.11). Substitution of (4.27) and (3.35) in equa- 
tions (4.23), (4.24) yields the following heavy flavour coeflicient functions 

■)2 ,,2 n2 



^5'=' (i^. ^) = { - 1" ^- + • (4.28) 

+ + j5o,qC(2)P,™ . (4.29) 
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Notice that in the above expressions the /i^-dependence has completely dis- 
appeared due to the special choice of the couphng constant renormalization 
scheme. 

Since we have now all renormalization group coefficients at hand we can 
calculate the heavy quark coefficient functions Hj^ i and j in the asymptotic 
limit ^ for arbitrary z. The splitting functions Pi?\ PiP and the 



nonpole terms in the OME's a^f , a\f = q,g) were presented in section 

3. The coefficients c'H, C2} {k = 1,2; i = q,g) appearing in the light quark 
and gluon deep inelastic coefficient functions can be found in Appendix B of 
[|rT[] (see also Appendix B in |2^). The analytic expressions for eqs. (4.14), 
(4.15), (4. 20), (4.21), (4.28), (4.29) are too long to be presented here and one 
can find them in Appendix D. 

Before finishing this section we want to study the behaviour of the coeffi- 
cient functions Hk,i{z, Q'^/m^, /rn?) {k = 2, L;i = q, g) in the limit z — >• 0, 
or 7] 00. The behaviour of these functions for z = at arbitrary has 
been predicted in pH] based on the method of /c^-factorization. Using the 



notations in [O we have the following predictions from pB 



2:^0 ' ^ TTL'^ 



(4.30) 



hm i/g) (z, e, 4) = -levraT^e [Gt{v, + Gl{v, 



(4.31) 



with i = q,g and Cg = Cp, Cg = Ca and r], ^ are defined in (2.17). The 
functions Gk{ri,^), Gk{rj,^) {k = L,T) are given by eqs.(19)-(24) in [|1^ for 
arbitrary ^. In the limit ^ 00 they behave like 



hm Gz.(r/,0 = 7^^[41ne-|] 

C^oo bTT 4 3 



11 14 14 

hm Gt{v, = ^7[2 In' e + ^ - 4C(2) + -] 

?-»oo bTT 4 3 3 



1 1, 



hm Gl{v,0 = 7--A-M, 
i^oo bvr 4 



(4.32) 
(4.33) 
(4.34) 
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lim (5^(7^,0 = ^^[-41ne + 2]. 



(4.35) 



Substitution of the above equations into (4.30) and (4.31) yields the asymp- 
totic behaviour 



2:^0 ' 2; L 3 3 /^^ 9 



(4.36) 



hmi/g)(z,e,^) = -C,7> 



16 

T 



1 2 > /32 , III ZUO\ , 



m 



208 ^ 



16, 32,, , 80 



(4.37) 



which agrees with expressions (D.3)-(D.6) in the hmit z 0. 
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5 Results 



In this section we want to make a comparison between the exact heavy quark 
coefficient functions in [§] and the asymptotic ones derived in this paper. The 
exact coefficient functions which are either available in computer programs 
[0] or in tables ]12|] were defined in the latter references using the notation 
— 2, iv; z = q,q,g;l = 0,1) and . These twelve coefficient 

functions are related to the Hj:\ in (2.14), (2.15) and L^i^\ (2.16) defined in 
this paper as follows. For the Born reaction we have the identity 

<i(^,0 = --cS(r],0, (5.1) 

with rj and ^ defined in (2.17). For the gluon-fusion process in (2.12) and the 
Bethe-Heitler process (2.13) we have 

= 167r^[cg(r^,0 +cS(r],01n4]' (5-2) 
and for the Compton process (2.13) one has 

4'K^'^'^) = 16vr^4>0. (5.3) 

The exact coefficient functions were plotted in figures 6-11 in as functions 
of 7] for various values of ^. Here we want to show at which ^ values the 
asymptotic forms of the coefficient functions presented in section 4 give a 
good approximation of the ones derived in 0. 

We do not need to discuss the lowest order coefficient functions c^?\ and 
C2°g as they have a very simple analytic form. The same is true for the coeffi- 
cients of the mass factorization logarithms, c^^^^^, C2^], cf^^ and C2^q. Therefore 
we concentrate on the asymptotic values of the NLO coefficient functions and 
start by defining the ratios b!]1\ and T^} which are given by 

<](^,e) = ^J^^, (5.4) 



and 



Tl^Kz,i)= fS^^'^^ . (5.5) 
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Here c'kJ'^^^'^i^^O ^"^^ d^kj^'^'^^^{z,C,) are the asymptotic expressions for the 
heavy quark coefficient functions derived in section 4 in the hmit ^ — >■ cxd. 
In this paper the above ratios will only be presented for charm production 
at the HERA collider on account of the small value of the charm mass and 
the values accessible. Choosing the range 5 < ^ < 10^ we will study the 
above ratios for z = 10~^ and 10""^. 

In fig. 5 we show R^^^g (5.4) and it is apparent that the asymptotic formula 
coincides with the exact NLO result when ^ > 10'^. There is essentially no 
difference between the ratios for z = 10~^ and 10~^. However if z = 10~^ the 
exact calculation begins to show computer instabilities when ,^ > 10^ so we 
did not run at larger values. 

The next figure fig. 6 shows R^^g and we deduce that our approximate 
formula is good for ^ > 10. This value is much lower than the one observed 
for R^^^g which shows that the asymptotic limit is reached much faster for ifg^j 

than for Hj^ g. Furthermore the lower bound is independent of our choice of 
z. The scatter in fig. 6 at large values of ^ reflects the numerical errors 



in the computer program of the exact NLO result in |T2|. This numerical 
uncertainty is therefore one of the reasons why we derived the asymptotic 
formulae for the heavy quark coefficient functions in this paper. 

Continuing with these ratios we plot R^^'^g versus ^ in fig. 7. There is not 

much difference between Ri^^g and R^^^g (see fig. 5) . The asymptotic formula 
for the quark channel is good for ^ > 10^ for both values of z. A simular 
observation holds for R'^]^ in Fig. 8 when compared with R'^g in fig. 6. Like 
in the latter case the exact formula approaches the exact one when ^ > 10. 
Finally we turn to the ratios T^^g (5.5) in Fig. 9. Here we have used 

(2) 

the exact NLO formulae for L), (5.3) presented in Appendix A rather than 



computing them from the program in |12|. This is the reason we have no 



numerical troubles at very large ^. In fig. 9 and fig. 10 we have plotted 
T^^g and Tg^^J respectively. From these figures we infer that at small z ( here 
z = 10~^ ) the asymptotic formulae coincide with the exact ones over the 
whole ^-range. At larger z-values (e.g. ai z = 10^^) the approximation gets 
worse and it becomes only good when ^ > 10^. 

Before drawing any conclusions about the validity of our asymptotic ex- 
pressions one has to bear the following in mind. First of all the heavy quark 
coefficient functions have to be convoluted with the parton densities accord- 
ing to (2.8) in order to obtain the charm contributions to the deep inelastic 
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structure functions. Therefore there are parts of the integration region over 
z which can become more important than some other ones. Second the lower 
bounds on the ^-values given above have to be viewed in the mathemati- 
cal rather than in the physical sense. Experimentally the structure function 
^2(3;, Q^, m^) will become much better determined than the longitudinal one 
given by Fl{x , , mf.) . In the latter case we are already glad that it can 
be measured up to 20 — 30% accuracy. Furthermore in this paper we are 
deahng with NLO corrections which have to be less precise than the Born 
approximation. Also the contribution coming from the Compton subproccss 
leading to the plots in figs. 9,10 is much smaller than the one originating 
from the photon-gluon fusion reaction for which the ratios are plotted 
in figs. 5,6. Therefore the physical bounds on ^ can be put much smaller 
than the ones given above. In view of the theoretical and experimental un- 
certainties we can state that the asymptotic formulae for the heavy flavour 
coefficient functions can be used when ^ > 4 in the case of F2 whereas ^ > 10 
is good enough for Fl- 

To summarize the calculations presented in this paper we have used the 
OPE techniques and the renormalization group equations to find analytic for- 
mulae for the asymptotic behaviour S> 1) of the heavy flavour coefficient 
functions which enter in deep inelastic electroproduction. We have tested 
these asymptotic formulae against the exact NLO results in our rather com- 
plicated computer programs and find good agreement when > 4 for F2 and 
^ > 10 in the case of Fl. Below these values our asymptotic formulae fail 
and one has to use the exact NLO computer programs to compute the heavy 
quark coefficient functions. 
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Appendix A 



In this Appendix we present the exact expressions for the heavy quark co- 
efficient functions L- ^ (2.16) corresponding to the Compton reaction (2.13). 
The calculation is straightforward because one can ffist integrate over the 
heavy quark momenta in the final state without affecting the momentum of 
the remaining light quark (see figs. 5c, d in 0). The phase space integrals 
are then the same as the ones obtained for the process 7* + g — > (7* + g 
{g* Q + Q) where the gluon g* becomes virtual. After integration over 
the virtual mass of the gluon one gets the expressions 

L^^liz, Q\ m") = CpTf [96z(|)2{ In + L.L^ + 2{-DIL, 

+DIL2 + DILs - DIU)] + )^(64 - 288z + lQ2z'^)Li 

rie 416/2X 1408 /2x2.^3 rl6 400 / 



3 V^y 3 Jsg2 ^3 18 ^{l-z)i 

160 3824 992 
9 



xf— + — -2-— 22)}sgiJ , (1.1) 



Lg(., Q\m-) = C,T,[{l'-±^ - ^(^'(1 -9z + 9z^)} 



X 



In 



Li + L1L2 + 2{-DILi + DIL2 + DIL3 - DILi)} 



3^ 1- z 



88 136 



(1 - 
152 1 



'(l28 -432z + 288^2 



9 



9 



-z — 



9 1-z V(l -2)^/V 9 



1 - z 

/464 512 2048 
-z + 



9 



832 6208 
+ —z - 



11392 2 , 6016 L3 



9 



9 



9 



-z^ + 



272 1244 718 1 ^ z 
~ ~rf ~ ^ 27 1-z ^ HI -2)^ 

4304 . 20 1 
— ^ + 27 



-z 



sq2 

3424 15608 
+ —r^Z 



27 



27 



[1.2) 
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where ^ = jrv? (2.17). Further we have defined 



sgi = Jl-4§ , sg2 = Jl-4 ^ 



L, = ln(l±^) , L. = lr.(l±^) , L3 = In (^^^) 



with < ^ < ^/(^ + 4). 
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Appendix B 



Here we apply the method of integration by parts which enables us to reduce 
scalar two-loop OME integrals with five diff'erent propagators to expressions 
containing only four different propagators. The method is a generalization 



of the trick invented in p3| where it was used to reduce two-loop self-energy 
integrals with five different propagators to integrals which only contain four 
different propagators. 

Feynman integrals with five different propagators emerge from the com- 
putation of the graphs e, /, h, I, m, n, o in fig. 2. Let us first start with graphs 
21, 2o which lead to the following expression 



j(m) 



(A ■ qy 



where a is 1 or 2. The denominators are defined as 



D. = k^ 



m 



1 5 



(2.1) 



(2.2) 



Dq = — m\ 



2 1 



D 



kp 



{k — pY — m\ 



Dqp = {q-p) - m: 



2 

■4 5 



Dkg = {k- q) - . 
In n-dimensional regularization the following integral is zero 







dJ'k d 



K{A ■ qy 



DkD-DkpDqpDkg 



(27r)"y (27r)"(9fc^ 
Differentiating the right-hand-side of (B.7) we obtain 

f d^k (A-g) 







{ 



(27r)"7 {271^ DkD-DkpDqpDkg 
2P 2k-{k-p) 2k-{k-q) 

Dk Dkp Dkq 



■}• 



(2.3) 
(2.4) 
(2.5) 
(2.6) 

(2.7) 



(2.8) 
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In the case that the following conditions are satisfied i.e., mi = 0, — 

ml + m|, ml = mf + m^ one can write 

-2k-{k-p)^-Dk-Dkp, (2.9) 



-2k-{k-q)^ -Dk - Dkg + Dg . 
Substitution of (B.9) and (B.IO) into (B.8) yields 



(2.10) 



1 



n 



(27r)"i (27r)"^ ^' ^D^Dl 



'^qK^^pDkq 
}. (2.11) 



If a = 1 we have already reached our goal. When a = 2 we have to repeat 
the trick for the last term in (B.ll) so that we finally get 



1 



n-AJ {2'nY J {2'kY 



(A ■ qy 



1 



D^qDlpD,pDkq 



1 



D-DkpD,,Dl^ 



n 



-{- 



5^D^q-'Dlp^,pDl, 



D--WkpD,,Dl^ DkD--WkpD,,Dl^ 



}]■ 



(2.12) 



For graphs 2e, 2h we have to differentiate as follows 

d'^q f dJ'k d r{k-q)^{A-qy 







J (27r)" J (2^ 



(27r)"i (27r)- dk^WkD-DkpD,pDkg 



dJ'q drk 



(A . qy 



X 



n/(™) + J 

r 2k-{k-q) 2{k-q)-{k-p) 2{k - qf 



(27r)" (27r)" DkD^DkpDgpDkq 



D 



kp 



D 



kq 



}. (2.13) 



Imposing the following conditions m^ — 0, m| 
one gets identity (B.IO) and 



ml + ml, ml = m\ + ml 



-2{k-q)-{k-p) = Dgp-Dkq-D_ 



'kp 



(2.14) 
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so that J^*"^ equals 

n-AJ 
1 

4 



{2t:Y J (27r)" 
1 



(A ■ q)' 



DlD-DkpD,, 
1 



DkD-Dl^Dk, 



DlD--WkpDgpDkg 



(2.15) 



For a = 1 we have aheady four different propagators. If a = 2 one has to 
repeat the trick again analogous to (B.ll). 

The trick of integration by parts also applies to graph 2m. For this graph 
we have the integral 



j(m) 



(A ■ qy 



d''q 

(2^ J (27r)- DkD-DkpDkgDkgp ' 



with 



D 



kqp 



{k-q-pf 



(2.16) 



(2.17) 



Performing the differentiation with respect to k in the same way as in (B.7) 
one obtains 



j(m) 



n 



i_ 

-47 (27r)"7 (27r)"^ ^' 

1 1 

+ 



+ 



D-Dl^DkgDkpg D-DkpDl^Dkgp 



1 1 
D.Dr'DkMLD 



(2.18) 



In the above expression we have used the relations in (B.9), (B.IO) and 



2k-{k-q-p) = 



D 



kq 



D 



kp I 



(2.19) 



2 2 n 2 2 

mg , m| = U, m2 — mi. 



which only holds under the conditions 

The integrals corresponding to graphs 2/ and 2n need some special treat- 
ment because of the sum X]i^o(^' ^)'(^'?)™ * which appears in the operator 
vertex. Notice that this sum also appears in graphs 2/i, 2o but after some 
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rearrangement of terms it drops out of the integrals. For graph 2n we have 
the expression 



d^'q f d''k (A ■ A;)™ - (A ■ q)' 



1 



(27r)-y (27r)- (A-k-A-q) D^DkDkpDkgD 



kqp 



, (2.20) 



with ml — ml = ml = m? and mi — mi — 0. The integral K^'^^ can be spht 
into X(-) = - with 



(27r)" 7 (27r)" (A • A; - A • g) D,DkDkpDk,Dk,p ' 



(2.21) 



(m) 



r d^'k (A ■ g)"^ 1 
(2^7 (27r)- (A • - A • g) DgD^DkpDkgDkgp 



(2.22) 



Differentiating the integrand in (B.21) with respect to q and (B.22) with 
respect to k analogous to (B.13) we obtain 



(m) 



1 



(A ■ ky 



n-bJ {2ttY J (27r)" (A • A; - A • g) 
1 1 



DkDgDkpDkgp 



DkDqDkpDj^^p D'^DkpDkqDkpq 



(2.23) 



(A ■ qy 



1 



n-Sy (27r)«y (27r)"(A-A;-A-g)LD2^^DfepL»fc,p 
1 1 1 



DqDlpDkqDkqp DkDqDl^Dkq DkDqDkpDl^p 
1 1 



DkDkpDkqDkqp DkDf^pDkqDkqp 



(2.24) 



As = the last two terms in (B.24) are equal to zero. 

Finally we apply the method of partial integration to graph 2/. The 
corresponding integral becomes 



d'^k (A • A;)™ - (A ■ qy 



1 



d^'q 

(2^7 (27r)'* (A-A;-A-g) DkDqDkpDqpDkq' 



(2.25) 
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with ml 



rrio 



ml. 



m^ and m^ = 0. L*^™) can be split in the same 



way as done for in (B.20). However because of the symmetry in k ^ q 
one can simphfy the calculation. Here we have 



^ J {2nY J 



(A ■ k)' 



1 



(27r)' 



-. (2.26) 



(27r)" (A ■ /c - A ■ g) DkDgDkpDgpDk, 
Differentiating the integrand with respect to q analogous to (B.13) we obtain 

1 



^(^) ^ _J_ r r d^k (A • kr 

^ n-bJ {2t:Y J (2tt> (A ■ k - A ■ 



(27r)" J (27r)" (A-k-A-q) iDkD^DkpDgp 
111 



DkDgDkpD^ D^^DkpDgpDk, D^D^Dl^D^^ 



(2.27) 
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Appendix C 

^ (2) 

Here we present the unrenormalized operator matrix elements A^j whose 
general structure expressed in renormalization group coefficients was derived 
in section 3. After having carried out mass renormalization the two-loop 
OME in fig. 2 is given by the following expression ( see also (3.22) ) 



A 



4{CfT/[(32 - 64z + 64^2) ln(l 



-(16 - 32^ + 64^2-) In ^ _ ^8 _ 32^)] 

+CATf [ - (32 - 6Az + Uz^) ln(l - z) - (32 + 128z) In ^ - 



64 
3i 



-16 



496 

128;^+— 



+■^{(7^7/ [(8 - 16z + 16z^)[2 lnzln(l - z) - W(l - z) + 2C(2)] 

-(4 - 8z + 16^^) In^ z - 32z(l - z) ln(l - z) 
-(12 - IQz + 32^2) In 2 - 56 + 116z - SOz^ 
+C^Tj [(16 + 32z + 32^2) [Li2(-z) + In z ln(l + z) 
+(8 - 16^ + 16^2) ln2(l - ^) + (8 + 16^) W z 



352 



+32zC(2) + 32^(1 - z) ln(l - ^) - (8 + 64z + ^ 



') In; 



160 
"97 



+ 16 - 200^ + 



H=c,b,t 

x(l4 



9 



-C(2) 



e/2 r_l 
e 



1744 (2) ■ 

rl 2/ 64 128 128 



(3.1) 



with 



j{2). 



;2) = CFTf[{l -2z + 2z^) [8C(3) + 8C(2) ln(l - z) + - ln^(l - z) 
-81n(l - ^)Li2(l - -z) + 4C(2) In^ - 41n2ln2(l - z) 

8 In ^Li2(l -z) + 8Li3(l - z) - 24Si,2(l - z)] 



2n 3 
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r 4 

- 24C(2) lnz + -ln^z+ 16 In 2^2(1 - z) + 32Si,2(l - z) 

-(4 + 96z - 64^2)Li2(l -z)-{6- 56z + 40z^)C(2) 

-(8 + A8z - 2Az^) In z ln(l - + (4 + 8z - I2z^) ln^(l - z) 

-(1 + I2z - 20^2) In^ z - (52z - 48z^) ln(l - z) 

-(16 + 18^ + 48z^) In z + 26 - 82z + 80^^} 

+CATf{{l -2z + 2;22)[_8((2) ln(l -z)-^ - z) 

+8 ln(l - 2)Li2(l - z)- 8Li3(l - z)] + (1 + 2^ + 2z^) 

X [-8C(2) ln(l + 2) - 16 ln(l + z)U2{-z) - 8 In z ln2(l + z) 

+41n2 ^ln(l + z) + 81n^Li2(-^) - 8Li3(-^) - 16Si,2(-;2)] 

4 8 

+(16 + 64^)[2Si,2(l -^) +lnzLi2(l -2)] - (3 + 3-^) In^-^ 

+ (8 - 32z + 16z^)C(3) - (24 + 96z)C(2) Inz + (16 + 16z^) 

^32 _ „, 272 
^3z 

16 



x[U2{-z) + ln;2ln(l + z)] + (— + 12 + 64;^ - —z')U2{l - z) 



-(l6 + 80z- 128z^ + — )C(2) -4^^ ln^ln(l-^) 

-(2 + 8^ - 10^2) W(l - ^) + (2 + 8^ + yz^) In^ z 

/. 2M /. N / 56 172 1600 2\, 

+(4 + 16z - 16^2) ln(l - 2) - (y + —z + ) In ^ 

448 4 628 6352 ^^ 
The unrenormalized OME corresponding to fig. 3 ( see (3.29) ) is given by 

64 ox 1/ , ,,2 / 64 px, 160 

+16^ + — / j + - (8(1 + z) W ^ - (8 + 40^ + —z^\ In ^ - — 

+16-48.+ l^.V«?/'^(-)}' (3-3) 

with 

a^f^^z) = Ul + ^)(32Si,2(l -z) + 161nzLi2(l - z) - 24C(2) In^ 
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+ ( _ 1^ _ 12 + 12^ + 16z^)c{2) + (2 + Wz + —z^) In" z 



z 

'56 88 



/uu oo 448 9\ , 
-{- + -z + —z^)lnz- 



448 



124 



1600 



27z 3 3 27 



3 3' 9 

The unrenormalized OME corresponding to fig. 4 is (see (3.33)) 

,2 ^2 



(3.4) 



iNS,(2)/A^ \ 2/"^ r 1 [32/ 1 \ 16 16 ■ 



1 

+- 

e 

+a 



9 VI 



8 1 + 8 88 

H m ^ H z 

3 1-z 9 9 



+ 5(l-.)(^C(2) + 9 



NS,(2), 



3> 

(3.5) 



4 4 N,,^, 1 + ^2 2 20, N 



with 

8, 224 / 1 N 44 268 

+'5(l-.-)(-5C(3) + fc(2) + g). (3.6) 
Here the term 1/(1 — z)+ has to be defined in the distribution sense 

/' dz{-^) f{z) = f dz-^[f{z) - /(I)] . (3.7) 
JO ^1 — z^ + Jo 1 — z 

Notice that as long as z < Q^/ (Q'^+Am^) with 7^ the terms proportional 
to 5(1 — z) do not contribute and the subscript + in (1/(1 — z))+ can be 
dropped. We will comment on this in Appendix D after having predicted 
the asymptotic expression for L^^g^^\z,Q^,m^). The OME's which emerge 
from the calculation of graphs 1 — 4 usually appear in the Mellin transformed 
representation 



A^;) = ^[i + (-i)-] I'dz z--'Aj{z) 



(3.8) 



This implies that the anomalous dimensions 7}™'* correspond to the physical 
operators listed in eqs.(2.24)-(2.26) only for even m. 
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Appendix D 

In this appendix we present the heavy quark coefficient functions HIJ and 
^f] — L,2;j — q,g) in the asymptotic hmit » m^. Starting with 
process (2.11) the heavy quark coefficient functions read ( see (4.6), (4.7) ) 

H^^l{z,Q\m')^Tf[16z{l-z)], (4.1) 



H^]l{z, g^ m^) = 7>[(4 -8z + 8^^) (in ^ + ln(l -z)-\n z) 

-4 + 32^ - 32z^] . (4.2) 

In next-to-leading order the coefficient functions given by process (2.12) and 
expressions (4.14), (4.15) read 

HQ{z,Q\m'^) = [CFT/{32^1n^ + 16 + 16^-32^2} 
+CAr/{64^(l - z) ln(l -z)- 128z In ^ + - 32 



544 

-160z + —z^} 



3z 

544 g2 

32 



In-^ 



+CATf{64:z{l - z) ln(l - z) - 128z In z + — - 32 

544 9^ , 

-imz + ^2 [ In — 

3 J 

+C^T,{(^ - + ipz3)[Li2(-z) +lnzln(l + z)] 

+32z[Li2(l -z)+\uz ln(l - z)] - (y z - ^^^)C(2) 

-f — z + —z^) In^ 2 + (16 + 48^ - 64^2) ln(l - 2) 
^3 5 ^ 

/ 64 208 416 192 2\ , 64 256 608 672 

+ (-15i-l^-— ")^"^+15^-l^-^^ + ^"l 
+CAT/{(64^ + 64z2)[Li2(-z) +ln^ln(l + 2)] - 128^112(1 - 2) 

+64^2^(2) - (192^ - 64^^) ln2ln(l - z) + {32z - 32z^) ln^(l - z) 

,39 928 \ 

+96^ In^ ^ + ( 32 - 288z + ln(l - z) 

^iiZ 3 ' 

+(32 + 256. - 416.^) In . - ^ + f + - l^.H , (4.3) 
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CFTf{{8 - 16z + 16^2) ln(l - z) 



-(4 - 8^ + 16^2) In ^ - 2 + 8^} 

+CATf{{8 - 16z + 16z^) ln(l - ^) + (8 + 32z) hiz + 



16 

3z 



124 

+4 + 32z - —z^} 



+ 



CfT/{(8 - 16z)Li2(l -z)- (32 - 6Az + 64z^)C(2) 

-(24 - A8z + 64z^) In z ln(l - z) + (16 - 32z + 32z'^) \n^(l - z) 

+(8 - 16z + 32z^) In^ z - (28 - 96z + 80z^) ln(l - z) 

+(8 - 48z + 80z^) In z + 36 - 68z + 16^^} 

+CAr/{-(16 + 32^ + 32^2) [U2{-z) +l-az ln(l + z)] 

+(16 + 64^)Li2(l -z)- (16 + 32z^)C(2) 

+(96z - 32^2) In z ln(l - ^) + (8 - IGz + 16z^) W(l - z) 



(16 + 48z) In' z + - 8 + 160z - ^z"") ln(l - z) 



536 



-(192z- 200z^)lnz + 



3;^ 3 

208 220 368 1628 o,i, 



9z 



-z + 



In 



32 



+CATf {(16 - 32z + 32^^) ln(l - z) + (16 + 64z) In^ + — + 8 
L 32; 

248 
3 

+(16 + 64^)Li2(l -z)- (16 - 32z + 32z2)C(2) 

+(96z - 32z^) ln^ln(l - ^) + (16 - 32z + 32^'^ ln2(l - z) 



+ 32z) In'z+i— 



3z 



lQ2z-^z^) \n{l-z) 



I 248 
-(8 + 256Z- — 

+Ci.r^ [l6(l + ( - 4Si,2(-^) - 4 ln(l + ^)Li2(-^) 

-2C(2) ln(l + - 2 In z ln'(l + + In' z ln(l + ^) - -Li2(l - z)) 

+8(1 - ^)'( - 4Si,2(l -z)- 3Li3(l - + 12Li3(-z) + ln3(l - z) 

In ^ ln'(l - ^) + 2 In' z ln(l - ^) + ln(l - ^)Li2(l - z) 



16 172 968 



z I In^ + —z + 



1124 



2 

, m 
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-41n2Li2(-^) +4C(2)ln^ - ^In^^ + 31n^ln(l - ^)) 
+8^'(4Si,2(l -z)- Li3(l -z)+ \n\l - z) 

7 

--lnzln2(l - z) + 41n2 ^ln(l - z) - ln(l - z)U2{l - z) 

-4C(2) ln(l - ^) + 41n zLi2(l - + 4C(2) Inz - In^ z + llLi2(l - ^)) 

+128^Li3(-z) + (112 - 96z + 192^2)C(3) - (112^ - 144^^) In z ln(l - 

/ 16 128 192 , , , , , 

+ (^^ + 96+ — ^+ — /)[Li2(-z) +ln2;ln(l + 2;)] 

208 1Q2 
+(48 - -^z + 104^2 + — ^^)C(2) - (22 - 88z + 84^^) ln'(l - z) 
3 5 

/ 8 96 \ 

-(4 - + 52^2 + — / ) In^ 2 + (28 - 132^ + 96^^) ln(l - z) 



3 5 

16 712 136 672 o\ , 16 904 68 328 



-I \ ln2;H \ z + 



15z 15 5 5 ^ 15z 15 5 5 

2n/^./1 — 1 — -^^ 



+C^T^ 16(1 + 2^ + 2^^) Li3(^) - Li3(-:^) 

'V 1 + 2; l + Z 

-ln2;ln(l - z) ln(l + ^) - ln(l - ^)Li2(-^) 

3 3 
+|ln^zln(l + ^) + -lnzLi2(-^)) 

+8(1 + 2;^) [5Si,2(l -z)+ 2Si,2(-;2) - 3Li3(-^) + In z W(l + ;2) 

+2 ln(l + z)Li2(-^) + C(2) ln(l + z)\ 

+32z[2Si,2(l - ^) + In ^Li2(l - z)\ 

-(16 + 128^)Li3(l - ^) + 8^^[ln^ zln(l + z) 

-2 In z\A2{-z)\ + (48^ - 16^^) In z ln^(l - z) 

-(8 + 64^ - 16z^) In^ z ln(l - z) + (16 + 64z) ln(l - z)Li2(l - z) 

-(40 - 48z + 64z2)C(2) ln(l - ^) + (y + 16z) ln=^ ^ 

-(16 + 160;^ - 32^2)C(2) Inz - (12 + 56z + 8z^)C(3) 

- (1^ + 48 - 16z - ^z^A [Li2(_^) + In ^ ln(l + z)\ + + 20 
^oz o ' ^oz 

-Uz + y^^)Li2(l -z)-i^-^ + A + 208z - ^-fz')a2) 

+(16 - 288^ + 292z^) Inz ln(l - ^) + - 6 + 64z - ^^^) ln2(l 

^oz 3 ' 
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.,,2 /208 112 860 2996 ox , , 
+(184^ - 114^2) W z + (— - —z + -g--2 ) ln(l - z) 

/292 5780 80 466 260 878 oi 

+ ( + 332^ z^]\^z^ — ^ + z z^ . 4.4 

V 3 9 y 9^ 9 9 9 J ^ ^ 

The coefficient function corresponding to the Bethe-Heitler process (2.13) 
reads (see (4.20), (4.21)) 

//^(z, g^ m^) = C^t\{ - 32^1n^ + I - 32 + ^^^j 

x{ln^ + ln— } -32^[Li2(l -^) + ln^ln(l -^) -In^^l 
+ (^ -32 + ^z^) ln(l -2) + (32 - 32^-64;2=^)ln;2 



^3^ 3 
32 32 128 320 

"9i + y-^"+^" 



(4.5) 



Q', m^) = CpTf [{8(1 + ;,) In ;^ + - + 4 - 4^ - -^2} In 



?>z 3 m 

+{[16(l + ^)ln^ + ^ + 8-8^-f ^2]!^^ 

+8(1 + ^)[2Li2(l - z) + 21nzln(l - z) - W z] 

+ (| + 8-8.-|.^)ln(l-.) 

/ 32 2\ , 16 160 16 128 2I , 

-(8 + 40.--.')ln.+ --- + -. + -.^}lu-^ 

+{16(l + z)[Li2(l -^) + ln2ln(l -z) -In^ z] 

32 32 2x, , X 2, 208 208 

+(— + 8 - 8^ - —z^) ln(l -z) + 32z^ \nz + 



2 



3z 3 ' ' ' 9z 3 

160 64 2. , 

H z 2;^}ln-+r 

3 9 m'^ 

+(1 + ^) (3281,2(1 -z)- 16Li3(l - z) + 81n^ W(l - z) 

-161n2zln(l -^) + 161n(l -z)Li2(l -z) -32C(2)ln^+ — In^ z) 
on QO 

-(— + 32 + 32^ + —z^) [U2{-z) +\nz ln(l + z)] 
^6z 3 ^ 

+ (^ + 16 - 16z + ^;.^)Li2(l -;.)-(- + 16 + 16z - ^z')a2) 
^oz o ^ z 3' 
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16 16 
+32^^ In z ln(l - ;z) + (— + 4 - 4z z^^ \^{\ 



+ (40z- 16z^)ln2z+ ( 



'iz 

208 208 



3 
160 



9z 



+ 



/280 



704 



- 88z z^] Inz 

3 9 ^ 

304 1216 832 



^9z ^ 9 



9 



-z + 



9 



(4.6) 



Finally we present the coefficient functions corresponding to the Compton 
process (2.13). They are given by eqs. (4.28), (4.29) which become equal to 



'NS,(2) 



z,Q\m^) = CFTf 
16 400 



}^z( ln^ + ln(l-2) -21n2 



Q 

m 



(4.7) 




gC(2)-yln^ln(l-^) 

4,2/, X .,2 58, ,^ , 134, 359 
+ - ln2(l -z)+A\n'z-— ln(l - 2) + — In ^ + — 

/2 26 X, ,^ , / 46 X, 29 295 1 

In the above expressions one should bear in mind that the singularity at 2; = 1 
will never show up because of the kinematical constraint z < Q'^/ (Q^ + 4m^). 
However after convoluting by the parton densities the structure function 
F2(z,(5^,m^) in (2.8) will diverge as ln^(Q^/m^) in the limit m^. 
In the above limit the upper boundary in (2.8) z^ax ^ 1 and the virtual 
gluon which decays into the heavy quark pair becomes soft. The soft gluon 
contribution which causes the cubic logarithm above has to be added to the 
two-loop virtual gluon corrections calculated in appendix A of . The 
cubic logarithm is then cancelled. The final result will be that in (D.8) the 
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singular terms at z = 1 have to be replaced by the distributions (ln*^(l — 
z)/{l — 2;))+ defined by 

and one has to add the following delta function contribution 



L2,q "{z,Q ,m)^CFTf5{l-z)[2ln - [^^^^ + T. 



12 2 

268 



+ 



9 



(4.10) 
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Figure Captions 



Fig. 1. a,b : One-loop graphs contributing to the OME Agg. 

c : Heavy quark {mjj > Q"^) loop contribution to A^^^ . The solid line 
indicates the heavy quark Q. 

Fig. 2. Two-loop graphs contributing to the OME A^^. The solid line in- 
dicates the heavy quark Q. Diagrams s and t with the external ghosts 
(dashed line) are needed to cancel the unphysical polarizations which 
appear in the sum of (3.43). Graphs u and v contain the external gluon 
self energy with the heavy quark {m% > Q'^) loop. 

PS ("2) 

Fig. 3. Two-loop graphs contributing to the OME Aq^'^ ' . The solid line 
represents the heavy quark Q whereas the dashed line stands for the 
light quark q. 

Fig. 4. Two-loop graphs contributing to the OME A^^'^'^\ The solid line 
represents the heavy quark Q whereas the dashed line stands for the 
light quark q. 

Fig. 5. R^l]y plotted function of ^ for fixed z — ^ (solid line) and 
z = 10~'^ (dashed line). 

Fig. 6. R2^l plotted as a function of ^ for fixed z — 10~^ (sohd line) and 
z — 10~^ (dashed line). 

Fig. 7. Ri^^q plotted as a function of ^ for fixed z — 10~^ (sohd fine) and 
z — 10~^ (dashed line). 

Fig. 8. R^2^g plotted as a function of ^ for fixed z = 10 ^ (solid line) and 
z = 10~^ (dashed line). 

Fig. 9. t£] plotted function of C, for fixed z — 10 ^ (solid line) and 
z = 10~'^ (dashed line). 

Fig. 10. T^l^ plotted as a function of ^ for fixed z — 10 (solid fine) and 
z — 10~'^ (dashed line). 
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Fig. 2 




Fig. 2 (continued) 
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Fig. 4 



